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ABSTRACT 

Aims. We investigate the origin of a flux increase found during a transit of TrES-1, observed with the HST. This feature in the HST 
light curve cannot be attributed to noise and is supposedly a dark area on the stellar surface of the host star eclipsed by TrES-1 during 
its transit. We investigate the likeliness of two possible hypothesis for its origin: A starspot or a second transiting planet. 
Methods. We made use of several transit observations of TrES-1 from space with the HST and from ground with the IAC-80 telescope. 
On the basis of these observations we did a statistical study of flux variations in each of the observed events, to investigate if similar 
flux increases are present in other parts of the data set. 

Results. The HST observation presents a single clear flux rise during a transit whereas the ground observations led to the detection 
of two such events but with low significance. In the case of having observed a starspot in the HST data, assuming a central impact 
between the spot and TrES-1, we would obtain a lower limit for the spot radius of 42000 km. For this radius the spot temperature 
would be 4690 K, 560 K lower then the stellar surface of 5250 K. For a putative second transiting planet we can set a lower limit for 
its radius at 0.37 Rj and for periods of less than 10.5 days, we can set an upper limit at 0.72 Kj. 

Conclusions. Assuming a conventional interpretation, then this HST observation constitutes the detection of a starspot. Alternatively, 
this flux rise might also be caused by an additional transiting planet. The true nature of the origin can be revealed if a wavelength 
dependency of the flux rise can be shown or discarded with a higher certainty. Additionally, the presence of a second planet can also 
be detected by radial velocity measurements. 

Key words, starspots - planetary systems 

1 . Introduction of the host star and compared it with other known planet-hosting 

stars. They found no chemical peculiarity compared to the other 

TrES-1 was the first planet transi ting a bright star that was dis- stars 
covered by a wide-field survey ( Alons oet al.ll2004T) . and it has 
become a frequently observed and well characterized exoplanet. 

It orbits with a period of 3.03 days a bright K0V star with an With the aim to acquire a very precise transit light curve 

apparent magnitude of V=11.79, approximately 1 57 pc away of TrES-1, our group observed in the year 2005 several transit 

from u s. A detailed light curve analysis done by IWinn et alj events with the Hubble Space Telescope (GO-10441). The most 

200 7]) allowed to l o wer th e errors of the star/planet parameters, relevant finding has been a flux increase during one transit, 

Charbonneau et all d2005l) obtained the first direct detection which is clearly a feature coming fr om the TrES-1 system itself , 

of light emitted from an exoplanet using the Spitzer Space for a preliminary description see ICharbonneau et all d2007l) . 

Telescope by an observati on of the TrES-1 secondary eclipse. Coverage of numerous further transits events was obtained be- 

ISozzetti et all d2004 120061) analyzed the chemical composition tween 2004 and 2007 with the IAC 80-cm telescope (IAC-80). 

The motivation for this coverage was for one a long-term study 

Send offprint requests to: mrabus@iac.es of transit-timing variations and later, a follow-up of the HST 

* Based on observations made with the NASA/ESA Hubble Space observations. The analysis of the flux-rises detected in these 

Telescope, obtained at the Space Telescope Science Institute, which is observations constitutes the main part of this paper, whereas the 

operated by the Association of Universities for Research in Astronomy, results from the timing study will be reported elsewhere. 
Inc., under NASA contract NAS 5-26555. 
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Both on TrES-1 (Charbonne au et ail 12007) and on other 
systems (H D209458, by ISilval 2003 bas ed on data from 
iDeeg et all 2001; HD189733 by lPont etal.1 2007), there have 
been reports of potential stellar spots detected from flux-rises 
during planetary transits. In starspots, intense magnetic fields 
can suppress the convection of heat to the surface, and hence, 
that area will radiate less light, causing a dark area on the stellar 
surface. Starspots c ould cover up to 55 % of the stellar surface 
dO'Neal et al.ll998l) and are generally 500 to 2000 Kelvin cooler 
then the surrounding area. The possibility to detect these during 
transits was first discussed by ISilval {2003). Their simulations 
showed that a starspot fully or partly occulted by the transiting 
planet leaves a short brightness increase in the light curve of a 
transit. The flux variation feature from a starspot depends on 
two parameters: Its duration depends on the transit path-length 
across the spot and across the stellar surface, on the rotational 
velocity of the parent star, and on the projected orbital speed 
of the transiting planet. The ratio of transit path-lengths may 
be interpreted in first order, assuming a central i mpact, as a n 
area ratio between spot and the transiting planet (ISilval 12003). 
The feature's amplitude for another, depends on the temperature 
difference between the normal stellar surface and the spot as 
well as on the fraction of the planet's shadow that grazes over 
the spot, unless the projected planet comes to lie fully within 
the starspot, in which case only the temperatures are relevant. 

For example, iPont et all (120071) found clear evidence of 
starspots on the star HD189733 by observing a transit of its 
planet with the HST at a high signal-to-noise ratio. They found 
an increase in flux of 0.001 mags and attributed this to the planet 
occulting a starspot. On TrES-1, there has been a previous 
ground -based search for spots during transits by IWinn et ail 
(2007) with negative results. They gave an upper 2 cr limit for 
spots of / < 4x10 ~ 4 , where / is the product of area ratio and the 
intensity contrast. 

In this paper, we present a second hypothesis for the origin 
of the short flux-rise observed with HST, occurring if two 
planets transit simultaneously, and one of the planets occults the 
other: In a system of two transiting planets with different orbits, 
the known one has a short period (planet A) and the hypothetical 
other one has a much longer period (planet B). Assuming that 
the transit of the short-period planet, i. e. planet A, is being 
observed and that the centers of both transits coincide more or 
less in time, then planet B can be considered as moving slowly 
across the stellar disc. When planet A is eclipsing planet B; 
it is causing in this moment a flux increase in the combined 
transit light curve. Due to the short observation span, the ingress 
and egress of the much longer transit events of planet B might 
neither have been observed nor detected, see Figure [6] for an 
arbitrary example light curve. For a light curve of two transiting 
planets, the flux variation depends on the area ratio and on the 
impact parameter between both planets, whereas the duration 
depends on the area ratio and on the respective orbital velocities 
of both planets. 

In Section 2 we describe the HST and IAC-80 observations 
and data reductions. We identified an increase in flux in one 
HST transit observation. Diagnostics to differentiate between 
the two origins of the flux rise are presented in Section 3, where 
we evaluate its color dependency from red to blue and then 
further investigate the flux variability of the other IAC-80 light 
curves. A discussion of the results will take place in Section 4. 



2. Observations and data reduction 

Table Q] gives an overview of the observations. The epochs 
refer to the planetary orbit; for the calculated mi d -transi t 
times we used the ephemeris from lAlonso et al.l (|2004), 
T c = 2,453, 186.8060+ 3.030065 x Epoch. Partial transits were 
observed in three 'visits' by the HST on Nov. 19 2004, Jan. 
19 2005 and March 29 2005. These observ ations followed the 
protocol established by Bro wn et"ail (|200l|) , where a complete 
transit light curve is obtained from phased overlay of data from 
several short 'visits' of the transiting system during slightly 
different phases. This phasing is necessary to overcome the 
interruptions due to Earth's occultation of the HST. For these 
HST observations the ACS/HRC instrument with the grism 
G800L was used. We collected spectra in the wavelength range 
between 5794 A and 10251 A. The exposure time was 95 s. 

Analysis of the HST images started with the pipeline- 
calibrated data, which subtracts zero point offset (based on 
over-scan mean), a two-dimensional bias and a dark (scaled by 
the exposure time) image. A wavelength dependent flat field, 
appropriate to position along the dispersion direction, has been 
derived by fitting a quadratic polynomial in wavelength to the 
eight broad-band flats over a range from 435 to 850 nm. A 
global sky level is computed by taking the median within a large 
region offset from the first order spectrum; this sky level was 
then subtracted from the images. A total of 20 bad pixels (out 
of an extraction box totaling 10700 pixels) were flagged as bad 
and linearly interpolated over, before proceeding. Cosmic rays 
are eliminated by sigma-clipping over stacks of 108 images on 
a per pixel basis, taking into account fluctuations of the flux, 
which are correlated with minor x-y image motions and changes 
of the focus. The pointing offsets in x and y are evaluated by 
taking first moments of the (nearly point source) zeroth order 
image, and focus changes are tracked by fitting Gaussians in 
cross-dispersion along the spectrum and averaging them. For 
photometry summing over the whole first order spectrum sums 
are taken out to an isophote, at which the per image counts (in a 
mean image) have fallen to 220 electrons (compared to a peak 
of 108000). Sums in 500 A bins use the same cross-dispersion 
limits, and contiguous blocks of columns most closely matching 
the desired ranges are formed. In addition a smaller bin about 
62 A wide is formed at H-alpha. 

Summation over the full spectrum results in total counts of 
about 4.5xl0 7 electrons per exposure. Normalization is done by 
taking the mean over all points in orbits 2 and 5 and dividing 
by this. This sum spans wavelengths from 5392 to 10650 A. 
Similar normalizations are applied separately to the wavelength 
restricted bin time series. De-correlations are performed by 
evaluating a linear least-squares fit to points outside the transit 
as a function of x-y offsets, the cross-dispersion spectrum 
width and the time to remove any linear trends in each visit. 
The resulting function is evaluated at all points including 
inside the transit and then divided out. Our final scatter in the 
time series of individual 500 A bands is increased by about 
9% over the sum of Poisson and readout noise. Any residual 
band-to-band correlated noise must be relatively small given the 
close approach to being Poisson noise limited. 

At eleven nights between 2004 and 2007, we also gathered 
rapid cadence photometry of TrES-1 transits with the IAC-80 
telescope at the Observatorio del Teide, Spain. All observations 
were done in the Johnson R filter. In the first two years, we 
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used a camera with a lk x lk CCD, whereas in 2006, it was 
replaced by a camera with a 2k x by 2k chip. The pixel scale of 
this new CCD is 0.305 ", resulting in a field-of-view of 10.25 
x 10.25 '. Out of the eleven IAC-80 transit observations, only 
the seven that cover full transit events were useful for further 
investigations. TableQ]shows details of the observations used. 

The images of each night were calibrated using standard 
IRAF pro cedures. We carried out aperture photometry with 
VAPH0T dDeeg & Dovld 1200 lb on the target and on five to 
eight comparison stars of similar brightness within the same 
CCD frame. In VAPHOT, differential photometry is obtained 
by dividing the target flux by the reference flux, which was 
constructed from a weighted ensemble of comparison stars. 
The flux was then converted into magnitudes and corrected 
for differential extinction by subtracting a parabolic fit to the 
off-transit parts of each individual light curve. We estimated 
the photometric error using the standard deviation outside the 
transit, see Table [TJ Figure [TJ shows the transit light curves of 
TrES-1 observed with the IAC-80. 

By visually inspecting the phased HST light curve (see 
Figure [2}, we can identify an increase in flux during the tran- 
sit on the night of November 19, 2004. The feature has a height 
of 2.7 mmags, whereas the standard deviation out of transit is 
0.18 mmags. Hence, this feature cannot be attributed to noise. 
The feature also is highly unlikely to be caused by instrumen- 
tal effects. No correlated variations in x-y pointing or focus are 
present at the time of the flux rise, i.e. nothing suggestive of an 
instrumental effect. During the several HST orbits, outside of the 
transit no features remotely similar to this have ever been seen. 
The short flux increase during the TrES-1 transit is ~ 10 min. 
long. This correspond to 1/15 of the whole transit of TrES-1. 
This kind of flux increase can occur when TrES-1, during the 
transit, passes over a darker area on the stellar surface. 

3. Analysis of the observed flux rise 

The two hypotheses presented previously, e. g. that the darker 
patch is caused by a) a cooler starspot or b) a second transiting 
planet would differ by these properties: 

- A starspot should show a flux-rise with a wavelength de- 
pendency, with a higher flux rise in the blue than in the red 
wavelength range. A flux rise from a mutual planet-planet 
occultation should not have any significant wavelength 
dependency. 

- A starspot is a temporary phenomena on time-scales of 
weeks to months; spots appear and disappear with different 
sizes at different positions. The planet-planet occultation 
would be a very rare event, except if the two planets are in a 
resonant orbit. 



We checked if there is any wavelength dependency of the 
flux rise feature during the transit of TrES-1. To this end, we 
created an HST transit light curve model by using the phased 
light curve of all three first order spectra HST observations and 
assigned a zero statistical weight to the data points of the feature, 
i. e. the flux rise was not considered for the following modelling. 
For trans it mod elling we used the FORTRAN routines from 
Gimenez (2006) and the simplex-downhill fitting algorithm 
dPress et al.l 119921) . minimizing x 2 - The free parameters were 
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Fig. 2. Phased light curve of all three HST visits using first order 
spectrum, each color represents a different day, i. e. red corre- 
sponds to November 19, 2004 , green to January 19, 2005 and 
blue to March 29, 2005. A flux increase is seen in the transit light 
curve on November 19, 2004. The black solid line is an example 
of a model for a simultaneous transit of TrES-1 and an arbitrary 
second transiting planet using UTM (see text Section|4|. In this 
example, the second planet has an arbitrary period of 300 days 
and radius of 4.25 R e - A nearly grazing impact between TrES-1 
and the second planet was simulated. 



the ratio between the stellar and planetary radii, k = R p /R t , the 
sum of the projected radii, rr = (R,/a) + {R p /a), and the limb 
darkening coefficient u\ and U2- Next, we used the color bins of 
the HST grism spectra and fitted the limb darkening coefficient 
in each wavelength bin, keeping rr and k constant, as obtained 
by the best fit of the first order spectra. The residuals are plotted 
in Figure [3] 

We used the model-subtracted light curves and estimated 
the slope of the flux increase in ea ch wavelength bin following 
the procedure bv lPont et alj (|2007), see Figure [3] The estimated 
slopes were then normalized with the slope estimated from the 
summed light curve. Figure [4] shows the estimated normalized 
slopes for each wavelength bin and their corresponding error 
bars. The slope in the H-a bin is a little higher; if real, this 
would indicate that TrES-1 occults a H-a rich region during 
transit. However, the light curve in H-a is also noisier and its 
error bars are higher, indicating that its higher slope might not 
be real. The slope of the bin with the central wavelength of 7750 
A is lower, but we attribute this to a single outlier. Also the bin 
centered at 10251 A has high error bars due to a noisier light 
curve. 

We used the above estimated slope values and fitted twice 
a linear polynomial of the form y — A + BA to them, where A 
is the central wavelength of the bin and y the normalized slope. 
In the first fit, we fitted for both parameters, A and B, and in the 
second one, we kept B at zero and fitted for A only. The first case 
indicates a wavelength dependency, the second case indicates a 
wavelength independence of the flux rise. For the first case we 
obtained a reduced^ 2 of 2.8 and for the latter case a reduced 
X 1 of 2.3. After applying the F-test, we found that there is no 
significant difference between both hypothesis, wavelength 
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Fig. 1. Corresponding to TableQ] IAC-80 observations of TrES-1 transits. The upper light curve is after subtracting a common transit 
model of TrES-1 . The lower curve shows the observed light curve with the solid (red) line being a boxcar smoothing and the arrows 
indicating possible low significance detection of similar flux rises as seen in the HST data. 
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Fig. 3. Phased light curve of all HST observations with the best 
fit transit model subtracted. From the bottom to the top are the 
light curves in the bins with central wavelength of H-a, 10251, 
9753, 9247, 8745, 8245, 7750, 7260, 6762, 6250, 5790 A. At 
the top is the sum over all wavelengths. The short inclined solid 
lines indicate the slope that was fitted to the flux rise. 

dependency and independence, within a 5% confidence interval. 

Furthermore, we used our IAC-80 data to check if these 
observations contain the same or similar features of flux-rises 
during transit. However, the IAC-80 observations are noisier 
and hence, similar features cannot be seen as easily as in the 
case of the HST light curve. Therefore, we created a best fit 
model for the transits from the IAC-80 observations, using the 
same approach as before, and subtracted the respective model. 
Then we used the flux rise event of the HST light curve as a 
template and applied a match filter detection against this event 
to the model subtracted IAC-80 observations, as follows: We 
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Fig. 4. Values of normalized fitted slopes for the flux rise in each 
wavelength bin and corresponding error bars. The dashed line 
is a two parameter linear fit and the dotted line a one parameter 
linear fit. The two parameter linear fit indicates a wavelength 
dependency and the one parameter fit a independence. 



shifted the template in time along the model subtracted IAC-80 
light curves, fitting in each shift for Mf,, = Must * f, where 
M HST is the magnitude of our template and / is a multiplication 
factor, which is also the fitting parameter. The cadence of the 
HST observations is lower than of the IAC-80 and, hence, we 
interpolated our template, using a spline interpolation. The best 
fit for / over time can be seen on the left side in Figure [5] and 
a histogram of / is shown on the right side. In the histograms, 
the value of / should concentrate around if there are no 
features, while a value close to 1 may mean a detection of a 
similar feature than observed with HST. Taking the off-transit 
standard deviation of / with ~ 0.27 into account, we consider as 
candidates for the detection of such a feature values of / > 0.8. 
Inspecting the histograms in Figure [5] we identified two low 
significant events on HJD-245000 = 3959.41 and 4271.60. 
These events are also indicated with arrows in Figures[TJand[5] 



4. Discussion and Conclusions 

We observed transits of TrES-1 in 2004 and 2005 with the HST 
and identified an intriguing flux rise during one transit. We 
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Fig. 5. Left side: Best fit / for each transit, where / is a multiplication factor of a model template, see text in Section [3] The solid 
(red) line below indicates the duration of the transit of TrES-1 and the arrows indicate a possible low-significance detection of a 
similar feature in the IAC-80 data. Right side: Histograms of /. The dotted lines give histograms for the out-of-transit part, the solid 
(blue) lines represent histograms of the transit part. 
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further followed-up transits of TrES-1 during the years 2004 to 
2007 from the ground, finding two further possible flux rises, 
and used the color information in the HST data to assess the 
origin of this feature. We can identify two possible causes, 
one is the possibility of a circular starspot occulted by TrES- 1 
during transit and the other is the detection of a second smaller 
transiting planet in an orbit larger than TrES-1. 

A good fit for the feature's shape can be obtained by using 
a circular dark patch on the stellar surface and neglecting 
deformation of a circular disk due to projection on a sphere. 
Assuming further that the dark patch and TrES-1 have a central 
impact on each other and using the orbital velocity for TrES-1 
of v = 140 km/s, then, in a time interval of 10 min., TrES-1 
crosses 84000 km, which gives us a lower limit for the diameter 
of the occulted patch. The real diameter depends on the stellar 
rotation (for the starspot hypothesis) or on the orbital velocities 
of the transiting planets (for the planet-planet hypothesis). 

The host star of TrES-1 is an active star dSozzetti et alJ 

2004) which is favouring the spot scenario. In that case, we 
calculated first the spot radius. Since the parent star of TrES-1 
is slowly rotating with a rotational velocity of Is - 1.08 km/s 
(Lau ghlin et al.1 120051) . corresponding to a rotation period of 
38.4 days, the stellar rotation can be neglected and under the 
assumption of a central impact between TrES-1 and spot, the 
spot radius is simply half the diameter as obtained previously, 
i. e. 42000 km. An effecti ve temperature of the stellar surface 
of 5250 K was obtained bv lSozzetti et al.1 (|2004) and combined 
with the equation 1 given by Silva (2003), we estimated the spot 
temperature to be approximately 4690 K for a starspot with the 
radius of 42000 km. These values are perfectly reasonable for a 
large sunspot. 

To generate a realistic model that shows that the observed 
HST transit light curve might as well arise from two occ ulting 
planets, we used "Universal Transit Modeller" (UTM, Deeg 
2008). UTM is an IDL-routine permitting simulations of transit 
light curves of multiple bodies, including rings and moons, and 
the fitting of observed light curves. We fitted first the summed 
HST light curve without considering the flux rise part. As fit 
parameters we used the radii of the host star and of the planet 
TrES-1, respectively; the limb darkening coefficient, assuming 
a quadratic limb darkening law, the inclination and the distance. 
We then used the best fit parameters and introduced a second 
transiting planet, but in an expanded orbit. However, we lack 
the determination of an unique period for the putative second 
transiting planet due to the singularity of the observed event and 
we further have no significant evidence in the IAC-80 data to 
derive a period. 

Hence, we can only give as an example a model of two 
transiting planets, see Figures [2] and [6] Figure [6] shows a 
complete example model, including ingress and egress of an 
arbitrary putative second planet. The introduced second planet 
has an arbitrary period of 300 days and radius of 4.25 R e . Its 
transit duration is 10 hours and the depth is 2.7 mmags. It is 
very unlikely that a complete transit of this second planet would 
have been observed, because the transit duration is too long for 
nightly observations and it is even more difficult to keep a suffi- 
cient high photometric precision from ground over such a long 
time span due to changing atmosphere, extinction and airmass. 
All IAC-80 observations were done in time spans shorter than 
10 hours and even under ideal conditions, i. e. at low airmass, 
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Fig. 6. Example of the same planet-planet transit model as shown 
in Figure [2] for TrES-1 and an arbitrary additional transiting 
planet, but showing the model over longer time span, so that the 
ingress and egress of a possible second planet are visible. The 
vertical solid lines indicate the usual observation span of 6 hours 
centered at the mid-transit of TrES-1. We note that during this 6 
hour time span the apparent transit depth of TrES-1 is identical 
with or without a second planet. 

the differential photometry applied to our transit observations 
did not correct completely for extinction. Therefore, as noted in 
Section [2] we corrected the extinction by subtracting a second 
order polynomial fit from the observed light curves which 
removed all possible hints to a second planet's transit from the 
light curve's off-transit part. 

Whereas a mutual occultation would be an extremely rare 
event, limits for the presence of transits from a second planet 
can be derived from size-estimates of this planet. Its minimum 
size can be obtained from the amplitude of the mutual transit 
event and assuming that TrES-1 and the other planet have a 
central impact on each other. In this case, the flux increase is 
directly related to the area of overlap of the projection of the 
two planets into the star, and is independent of the period. For 
TrES-1 we have a transit depth of 23 mmags and the flux rise 
height is 2.7 mmags, and hence we obtain an area of overlap 
of 1/8.5 of the cro ss section of Tr ES-1. Knowing the radius of 
TrES-1, 1.081 R y dWinn et al.ll2007l) . we obtained a minimum 
size of a second planet of 0.37 R/. However, the shape of the 
flux rise indicates a more grazing impact between TrES-1 and 
the hypothetical second planet, favouring a bigger radius. 

The original STARE and PSST data of the field where 
TrES-1 has been d iscovered cover a time span of 96 days 
dAlonso et al.l 120041) . This is a time span long enough to 
possibly detect a putative second planet on a longer orbit, but 
no transit has been detected in this data. Therefore, based on 
an assumption of a standard deviation of 10 mmags in the 
TrES-1 discovery data, we can set an upper limit for the size 
of a possible second planet. In the limiting case, a transit with 
a depth of 10 mmags would not have been detected. This 10 
mmags corresponds to a transit depth of a planet with a radius 
of 0.72 Ry, which constitutes an upper detection limit for the 
radius of a possible second planet. 
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For a putative long-period outer planet and due to the lim- 
ited and gaped observations of the TrES-1 field the detection 
depends mostly on the probability of having observed a transit. 
Therefore, we further constrain the period in which this upper 
radius limit is valid. We applied a statistical analysis to a time 
vector which contains the time of observation. We established 
2000 trial periods between a period of 4 days and 100 days 
on the STARE/PSST discovery data. Within each trial period 
we set randomly 300 epochs and tested how often we could 
recover a transit signal. We set a detection limit for a possible 
planet at periods where the recovery rate drops below 95 %. For 
the original STARE and PSST data this occurs at a period of 
10.5 days. Hence, for periods of less than 10.5 days, we can set 
an upper limit for the radius of 0.72 R j for a possible second 
planet. However, larger planets in longer periods are possible. 

Based on the conventional interpretation of flux-rise fea- 
tures during transits, we may conclude that a starspot has 
been detected in the HST data and possibly as well in the 
ground-based IAC-80 data. We also presented an alternative 
hypothesis of a mutual planet-planet occultation, which would 
be a rare but not impossible event. The verification or exclusion 
of this hypothesis is difficult unless a wavelength-dependency of 
the flux-rise features can be demonstrated, or until the presence 
of the secondary planet can either be verified or excluded, be 
it from observations of its own transits or from radial velocity 
measurements. 
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